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ABSTRACT 

We study possible astrophysical and dark matter (DM) explanations for the 
Fermi gamma-ray haze in the Milky Way halo. As representatives of various DM 
models, we consider DM particles annihilating into W~^W^, bb, and e~^e~. In 
the first two cases, the prompt gamma-ray emission from DM annihilations is 
significant or even dominant at > 10 GeV, while inverse Compton scattering 
(ICS) from annihilating DM products is insignificant. For the e~^e~ annihilation 
mode, we require a boost factor of order 100 to get significant contribution to the 
gamma-ray haze from ICS photons. Possible astrophysical sources of high en- 
ergy particles at high latitudes include type la supernovae (SNe) and millisecond 
pulsars (MSPs). Based on our current understanding of la SNe rates, they do 
not contribute significantly to gamma-ray flux in the halo of the Milky Way. As 
the MSP population in the stellar halo of the Milky Way is not well constrained, 
MSPs may be a viable source of gamma-rays at high latitudes provided that there 
are ~ (2 — 6) x 10^ of MSPs in the Milky Way stellar halo. In this case, pulsed 
gamma-ray emission from MSPs can contribute to gamma-rays around few GeV's 
while the ICS photons from MSP electrons and positrons may be significant at 
all energies in the gamma-ray haze. The plausibility of such a population of 
MSPs is discussed. Consistency with the Wilkinson Microwave Anisotropy Probe 
( WMAP) microwave haze requires that either a significant fraction of MSP spin- 
down energy is converted into e'^e~ flux or the DM annihilates predominantly 
into leptons with a boost factor of order 100. 
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Recently, iDobler et al.l ( l2010l ) have found evidence for a 7-ray haze in the halo around 

lis signal can be a signature of dark ma tter (DM) 



the Milky Way G alactic center (GC). 
annihilation (e.g. 



Zeldovich et al 



1980 



Springel et al 



I2OO8 



Kuhlen et al. 



2OO9I). The pri- 



mary purpose of this paper is to look for possible astrophysical sources of the haze and 
compare them with DM. Annihilating DM particles that produce many prompt 7's (via 
channels such as XX ~^ W^W~ , ZZ, bb, r+r~), may significantly contribute or even domi- 
nate at photon energies above 10 GeV. The DM particles that predominantly annihilate into 
leptons contribute significantly to the 7-ray haze only if their annihilation cross section is 



enhanced by a bo ost factor of order 



feld enhancement ( Hisano et al. 



2005 



00. This boost factor can be attributed to Sommer- 



Arkani-Hamed et al. 



annihilation channel yy — )■ (pep — )■ 2e^2e (iFinkbeiner &: Weiner 



20091) in, e.g . , XDM models with 



2007 



Cholis et al 



l2009b 



Arkani-Hamed et al 



20091) 



A useful discrimination of various sources is their total power in the Milky Way halo. 
In order to estimat e the power o f the 7 -ray haze, let us first calculate it in the 7-ray haze 
"window" used by iDobler et al.l (120101 ) to find the spectrum of 7-rays in the haze. This 
region is —15° < / < 15° and —30° < b < —10°. The corresponding solid an ^le is flh^vp 
ih. — M (sin&^ — sin 61) ^ 0.17. Integrating the 7-ray spectrum in Figure 11 of 



Dobler et al. 



mm . we find 



7— haze 



10^'^ erg s 



dE 



'V. 



where Rq = 8.5 kpc is the distance from the GC to the Sun. The haze is observed within 
approximately 9 = 45° from the GC. The corresponding solid angle Qtot = 2it{1—cos6) ~ 1.8. 
Therefore, the total luminosity of this 7-ray emission is 

fit 



'tot 



7 tot 



haze 



7— haze 



10^^ erg s"^ 



(2) 



At first we will discuss possible astrophysical sources of the 7-ray haze. The current 
star formation rate in the halo of the Milky Way is very small. Thus the sources of the high 
energy particles should have a long lifetime. The two possibilities are type la supernovae 
(SNe) and millisecond pulsars (MSPs). 

Let us estimate the output in the high energy electrons from the type la SNe. On 
average, SNe must produce ~ 10'^^ erg in relativistic electrons to account for the observed 



3 



flux of cosmic-ray electrons (IKobayashi et al.ll2004r). The ca . 



diet s imilar or smaller power in electrons ( iBerezhko fc Volk 



culat i ons for observed SNe pre- 



20081: 



Zirakashvili fc Aharonian 



20101). The birth rate of la SNe per unit stellar mass in the Galactic halo can be estimated 
as (5.3 ± 1.1) X 10~^^ yr~^ (ISuUivan et al.ll2006l ). In order to estimate the mass of the 



Mi 
by 



ky Way ste l 



Juric et al. 



ar ha lo, we use the distribution of matter in the disk and in the halo given 



( I2OO8I ) and, for the overall norrnalizat ion, we use the local stellar density of 
the thin disk, 35 Mq pc~^ (IKuijken fc Gilmord 119891 ). Therefore, the inferred mass of the 



Milky Way stellar halo within 20 kpc of the GC is 



Mhalo ~ 10' M( 







(3) 



with an uncertainty of at most a factor of 2. This gives a la SNe rate in the halo of about 
5 X 10~^ yr~^ or 2 x 10^^^ s^^. Consequently, the electron output of the halo la SNe is 



± la 



< 10^^ erg s"\ 



(4) 



This is insufficient to account for the 7-ray haze already based on the total output energy 
(we need at least 10^^ erg s~^). Provided that the above estimations are good within an order 
of magnitude, we conclude that la SNe will not contribute a significant number of electrons 
and 7-rays in the haze region and we will not consider them in the following. 



MSPs are known to emit pulsed 7-rays with a cutoff at a few GeV ( lAbdo et al.ll2009al Jbl). 
Observations o f X-ray nebula around some MSPs show a p roduction of high energy electrons 
and positrons 



( jstappers" 



et al. 



2003 



Kargaltsev et al.ll2006l ). but the uncertainties in the par- 



ticle spectrum are rather large. We will consider two possibilities for the total energy output 
in e^e~ . In the first case, we assume that the electron emission from MSPs is similar to the 



electron emission from young radio pulsars, such as the Crab pulsar ( lAtoyan &: Aharonian 



19961 ). In particular, we assume that a large fraction of the spin-down energy of an MSP goes 
into electrons and positrons with a cut-off in the e^e~ injection spectrum, E^nt ^ lOOGeV. 
In this case, we demonstrate that the MSPs may be sufficient to explain 7-ray haze at all 
energies (although prompt 7-ray emission from DM annihilation may contribute significantly 
at > 100 GeV). In the second case, we assume that e^e~ emission from MSPs is small, 
then the spectrum above 10 GeV requires an additional source, such as annihilating DM. 
At the m oment both possibilit ies for e^e~ emission are consistent with observations of MSP 
nebulae (IStappers et al.l 120031 ). The average prop erties of pulsed 7-r ay emission from MSPs 



for a sample of eight pulsars detected by Fermi ( lAbdo et al. 



2009al ) are presented in Table 



-4- 



Index r 


Cutoff Ecut (GeV) 


logio^ (erg s ^) 


1.5 ±0.4 


2.8 ± 1.9 


33.9 ±0.6 



Table 1: Some average properties of eight MSPs observed by Fermi (jAbdo et al.ll2009al l. Here we 
assume that the 7-ray spectrum is a power-iaw with an exponential cutoff, F ~ S~^e~^/^''"*. The 
total power in 7-rays equals the total spin-down luminosity times the conversion efficiency = rj^L. 
Assuming the log-normal distribution of L, the mean luminosity is (L) ~ 2 x 10^^ erg s~^. Assuming 
T]^ ~ 10% we get the mean luminosity in pulsed 7-rays to be (L^) ~ 2 x 10'^^ erg 



Let us estimate the number of active MSPs in the Milky Way halo. One of the standard 
models for the formation of the Galaxy assumes that the halo is formed before the disk 
(IBinney fc Tremaindl2008l ). The star formation in the halo happens predominantly during 
the early stages of the evolution. Based on the location o f observed MSPs on the PP 



diagram (ILorimer fc Kramer 



2005 



Chen fc Ruderman 



19931), the length of time they can 



produce e^e pairs in their magnetosphere is on the order of the Hubble time. Therefore, 
we expect that any MSPs produced during this period are still activ e today. This conclusion 



is also consistent with the o 



Ferrario fc Wickramasinghe 



Dserved characteristic age of the MSPs (IManchester et al. 



2007h . 



2005 



Since the direct detection of MSPs in the Milky Way halo is difficult, we need an indirect 
way to estimate their number. An important observation is that the number density of MSPs 
depends sensitively on the stellar density. For instance, globular clusters have much higher 
stellar density thari the galaxies. The 47 Tucanae cluster is expected to have about 50 MSPs 
( Abdo et al.l 2009b ) and a mass of 10^ M©, which gives the MSP number to stellar mass ratio 
of 5 X 10~^ M^^. For the Milky Way dis k, the local column density of MSPs is estimated 



to be 50 kpc fjCordes fc Chernoffl 119971 ). For a local total mass density of 50 Mq pc 



flKuijken fc Gilmorelll989f L ;his gives the MSP number to mass ratio of 10" 
50 times smaller than the corresponding ratio for the 47 Tuc. 



^ Mq^ which is 



The halo of the Milky Way is formed from both the initial s pherical stellar popu . 



and via mergers with nearby dwarf g alaxies and stellar clusters ( IBuUock fc Johnston 



ation 



2005 



Font et al. 



2006 



Zolotov et al. 



2OO9I ) . Galaxies at high redshift have significantly higher 
stellar density tha n the low redshift gala xies, for instance, sX z = 2 the density may be up to 



100 times higher ( Bezanson et al.l 2009). For comparison, the Milky Way stellar dei isity is 
of or der 10^ times smaller than the stellar density of 47 Tuc within half-mass radius (iHarris 
19961 ). We expect the MSP number to mass ratio to be determined by the stellar density at 
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the time of star formation. Thus the MSP number to stellar mass ratio for the Milky Way 
halo should be higher than that for the Milky Way disk but lower than that in the globular 
clusters. For a Milky Way stellar halo mass of Mhaio ~ 10^ M©, the expected number of 
MSPs is 



halo MSPs 



1 X lO'^ - 5 X 10^ 



(5) 



where the lower bound comes from the MSP number to stellar mass ratio for the Milky Way 
disk and the upper bound comes from the MSP number to stellar mass ratio for the 47 Tuc. 
Using the parameters from Table [H we find that the total power emitted by the halo MSPs 
can be as high as 

W^haio MSPs -10=^' erg s-^ (6) 

and the power in pulsed 7-rays can be 10^* ergs~^. This estimation can vary significantly 
due to uncertainty in the total number of MSPs in the stellar halo. It may also depend on 
the distribution of pulsar properties, e.g., the index, the cutoff and the luminosity in pulsed 
7-rays. With only eight observed 7-ray MSPs, this distribution cannot be well measured, 
but it should be possible in the future with a better observational statistic. Thus, at the 
moment, we cannot robustly estimate the contribution of MSPs to high latitude 7-rays. The 
main purpose of this work is to point out the possibility of a significant population of MSPs 
in the stellar halo of the Milky Way which are a viable source of the 7-ray haze. In this 
paper, we will take the 7-ray luminosity of MSPs in Table [1] as a reference value. Then, for 
some particular DM models, we will be able to estimate the number of MSPs in the stellar 
halo from the 7-ray haze data. 

MSPs are also continuously created in the Galactic disk. Every pulsar acquires a 
kick velocity at the birth. For r egular pulsars the mean kick velocity is about 400 km s~^ 



mass 



1997h . 



( iFaucher-Giguere fc Kaspil 120061 ). MSPs are usually found in bin aries. The center o : 
velocities o f the b inary systems have typical values ~ 80 km s~^ (jCordes fc Chernofl 
Lyne et al.l ( 1l998[ l obtained the average velocity ~ 130 km s^^ for their sample of MSPs. 
These velocities are not sufficient for t he majority of MSPs to overcome the gravitational po- 



tential of the disk (ICordes fc ChernofI 



1997 



Storv et al.ll2007h. Although the disk MSPs con- 



tribu te significantly to the diffuse 7-ray background ( IStorv et al. I l2007i iFaucher-Giguere Sz Loebl 
20101 ). the corresponding flux of 7 -rays has a disk-lik e morphology which is different from an 
egg-like shape of the 7-ray haze (IDobler et al.ll2010l ). Consequently, we need an additional 
population of MSPs that stretches to higher latitudes. 



The distribution of the stellar mass in the halo falls approximately as r ^ ( jJuric et al. 
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2008 



Bell et al. 



20081 ). In order to get a smooth distribution in the Galactic Center we 



introduce a cutoff radius -Re, and talce the following profile of MSPs 

1 



PMSPs 



ir + Rc 



(7) 



where we will take Rc = 2 kpc as our reference value. Physically, at R^ the power-law 
distribution of the matter in t he halo breaks do wn and merges to a Gaussian or exponen- 



tial distribution near the GC fiDwek et al. 



19951 ). Since we are interested in high latitude 



fluxes, the precise shape of the matter distribution near the GC is not important and the 
approximation in Equation ([7]) is sufficient. 

In our calculations, we will use the following source function for e"^ and 7 emission from 
the MSPs (the parametric form of the source functions for and for 7 are the same but 
the index and the cutoff are different) 



Qmsps('^ 5 E) — Qo- 



1 



r + Rc 



E 

Eriti 



exp 



E 



Er 



where for the 7-ray spectrum we will use the parameters from Table [T], n = 1.5 ± 0.4 
and Ecut = 2.8 ± 1.9. Motivated by possible similarities between the electron spectrum from 
MSPs and younger radio pulsars, we use the following parameters to de scribe the spectrum of 
electro ns: n = 1.5±0.5 and Ecut equal to a few hundred GeV (see, e.g.. lAtoyan &: Aharonian 
( 119961 ) for the e^e~ spectrum in the Crab Nebula). The total energy output is proportional 
to the total stellar mass of the halo in Equation ([3]) times the MSP number to stellar mass 
ratio in Equation (|5]). 



Annihilating DM provides a source for prompt 7-rays and e^e of the form 

QBM{r,E) 



1 Pdm dN 
2K^^''''^'dE-' 



(9) 



where Pdm/"^dm number density, (crf)o = 3.0 x lO^^^cm^s"^ is the thermally 

averaged annihilation cross section at freeze out, and dN/dE is the differential energy spec- 
trum of either or 7 produced in a single annihilation event. In general, the annihilation 
rate can be enhanced by a boost factor which we denote as BP. For the calculations we use 
the Einasto profile of DM dMerritt et alll2005 ): 



pDM[r) = PDMoexp 



2 



r: 



■0 



a 



(10) 
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where i?_ 



25 kpc, a = 0.17, and = 8.5 kpc. For the l ocal DM density, we take 



0.4GeVcm-3 Jcatena fc UlliJboid ). 



Merritt et al. 



( 20051 ) motivated the use of this 



PdMo 

profile by the DM A^-body simulations. The rotation curves in the inner part of the Galaxy 
are d ominated by the luminous matter a nd do not con s train significantly the DM distribution 



Dehnen &: Binneyi (119981 ) . see also 



Salucci et al. 



(120071 ) for a recent fit with a less cuspy 



(.e.g. 

profile). In any case, the spectrum of photons in the 7-ray haze window is not sensitive to 
the choice of DM profile. 

The total power released by annihilating DM is 



DM 



'IV. 



For a DM particle with the mass Mdm = 300 GeV and the rest of the parameters described 
above, the total power from annihilating DM within 20 kpc from GC is 



W^DM ~2 X 10^^ ergs-\ 



(12) 



which is about 5 times smaller than the total power in the 7-ray haze. Based on this crude 
estimate we already expect that without boost factors this DM model will not explain the 7- 
ray haze at all energies but it can be significant at some energies provided that the production 
of high energy photons is efficient. 



For the propagation of in the interstellar medium we use GALPRO P (IStrong et al. 



2007 



Moskalenko fc StrongI 1999 ). with the interstellar radiation field model of Porter fc Strong 



( I2OO5I ). We assume the following magnetic field model 

Rq — R 



B{R, z) = Bq exp 



Rr 



exp — 



(13) 



where Bo = 5 fiG is the value for the total (combined random and large scale "ordered") 
local magnetic field. Rc = 4.5 kpc and Zc = 2.0 kpc are the characteristic scales along R and 
z directions, respectively, and Rq = 8.5 kpc. 

Our main result is shown in Figure [H In order to explain the 7-ray haze, we add the di- 
rect emission of pulsed 7-rays from MSPs, the ICS photons from MSP and DM electrons, and 
prompt 7-ray emission from annihilating DM. The prompt 7-ray emission from annihilating 
DM becomes comparable to ICS photons from MSP electrons around 100 GeV. The bump 
in the 7-ray haze spectrum around 4 GeV is naturally explained with the pulsed 7-ray emis- 
sion from MSPs. We also compare the synchrotron radiation with the microwave haze data 
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Fermi Haze Template < 15, -30 < b < -10) 
Total 

DM W*W", Prompt j 

DM W*W\ ICS off ISRFxIO 

MSPs pulsed T, P,,=5.6x10^^erg s'^ 

MSPs ICS off ISRF, Pe==2.7x10^''erg s"' 



300 GeV 




Total 

DM W+W, M^=300 GeV 
Milky Way Halo MSPs 
constant 
Haze Data i 




10 15 20 25 30 

Latitudinal Radial Distance from GC (degrees) 

Fig. 1. — Left: contribution of pulsed 7-ray emission from MSPs, prompt 7-ray emission from 
annihilating DM, and ICS off MSP and DM electrons to the 7-rays haze spectrum. Right: cor- 
responding contribution of the synchrotron radiation from the electrons to the microwave haze at 
23 GHz. The parameters for the pulsed 7-ray emission are = 1.3, and Ecut^ = 4 GeV. The 
MSP e~^e~ injection spectrum parameters are Ue = 1.3, and i^cutc = 300 GeV. The DM has a 
mass Mdm = 300 GeV and annihilates into W^W~ with {av)o = 3.0 x 10"^^cm^s~^ The spatial 
distribution of MSPs and DM is discussed in the text. We use Rc = 2 kpc for the distribution of 
MSPs. The total power in pulsed 7-rays and in e~^e~ emission from MSPs is 



5.6 X 10^'' erg s ^ 



and W„ 



2.7 X 10^^ ergs"\ 



respectively. For a mean Lmsp ~ 2 x 10^^ erg s 



about 3 X 10'^ halo MSPs with average conversion efficiencies rjy ~ 0.1 and rj^ 



it corresponds to 
! 0.5. 



(jFinkbeineii |2004| ) . We add the synchrotron from MSP and DM electrons together with a 
constant a priori unknown offset. Our model seems to give synchrotron radiation that is less 
cuspy than the microwave haze data. The deviation happens below 10° which corresponds 
to a distance of about 1.5 kpc from the GC where our mass distribution in Equation ([7]) may 
be invalid. We also expect some contribution from sources in the bulge which may provide 
cuspier latitudinal profile of synchrotron radiation than the one presented in Figure [H 

In this example, we need about 3 x 10^ MSPs in the Milky Way halo which corresponds 
to (1 — 2) X 10^ MSPs in the 7-ray haze "window". The question is whether it is possible to 
distinguish between the 7-ray flux from MSPs and the diffuse ICS photons. The flux from 
an average MSP with pulsed 7-ray luminosity ~ 10'^^ erg s~^ at a distance of 8.5 kpc from 
the Earth is < 10^^^ erg s~^cm~^, which is less than the Fermi sensitivity for a detection of 
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Fermi Haze Template < 15, -30 < b < -10) 
Total 

DM W*W", Prompt j 

DM W*W\ ICS off ISRFxIO 

MSPs pulsed T, P,=1 .3x10^^erg s'^ 



M„=300 GeV, B.F.=3 




Total 

DM W+W, M^=300 GeV, B.F.=3 
constant 
Haze Data : 



V = 23 GHz 



10 15 20 25 30 

Latitudinal Radial Distance from GC (degrees) 



35 



Fig. 2. — Same as in Figure [U but witliout MSP electrons. For the 7-ray haze, we need about 
2.3 more power in pulsed 7-ray emission and a boost factor of 3 for the DM relative to Figure [TJ 
For a mean luminosity in pulsed 7-rays L^msp ~ 2 x 10'^^ ergs^^, it corresponds to a population 
of about 6 X 10^ MSPs in the Milky Way halo which is slightly higher than the upper bound in 
Equation ([5]). 

a point source with 80" significance even after a year of observationll] One can also use more 
sophisticated methods to distinguish between the point sources and the diffuse background 
without point source identification ( Slatyer fc Finkbeiner 20icl ). The problem is that the 
expected typical separation between MSPs in the 7-ray haz e "window" i s comparable to 
the Fermi-LAT point spread function at energies < 10 GeV (iRandd l2009l ). Consequently, 
it seems implausible that one will be able to distinguish between the 7-ray emission from 
MSPs and the diffuse ICS photons using available 7-ray observations. 

Let us argue now that a population of 3 x 10^ MSPs in the Milky Way halo is not 
inconsistent with the properties of currently known MSPs. Since the luminosity of MSPs 
is generally very low, most of the observed MSPs are within ~ 1 — 2 kpc from the Earth. 
Thus, we cannot directly observe the halo MSPs at large heights above the disk and the 
only parameter that can effectively distinguish a halo MSP from a disk MSP is the velocity 
relative to the Sun. Base d on the force in the vertical direction at R = Rq, z = 1.1 kpc 
( iKuijken fc Gilmord Il99ll ). K^^li = 2G x (71 ± 6)MqPc~^, we can estimate the escape 
velocity fesc ~ 100 kms~^. Pulsars with the transverse velocity Vtr ^ Vesc can be attributed 



^ Fermi-LAT performance: http: //www-glast . slac . Stanford. edu/software/IS/glast_lat_performcince .htm 
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to the halo population. 



According to the ATNF catalog ([Manchester et al.ll2005[ ). there are about three MSPs 
with vtr > 200 km at a distance r ~ 1 kpc from the Earth (e.g., B1957+20, J1909-3744, 
and B 1257+ 12). This is consistent with the local number density that follows from the 
distribution in Equation ([7]) nhaioMSP ~ 3 kpc^'^. One should also take into account that not 
all local MSPs can be observed due to beaming of the radio emission. Th e local MSP number 



Cordes fc Chernofi 



densit y that we need is significantly larger than the one estimated by 
(jl997l ). nh < 0.4kpc~^ at 90% confidence. We believe that this discrepanc y may be due to the 
small velocity dispersion = 53 km s~^ in the sample of MSPs used by 



Cordes &: Chernofi 



(119971 ). This velocity dispersion implies that the probability to have Vtr > 200 km s ^ is less 



than 10 ^ which is inconsistent with a few MS 



that have v^^ > 200. The analysis of iLyne et al 



^ s out of about 100 in the ATNF catalog 
^19981 ) gives = (80 ± 20) km s^^ which 



corresponds to the probability P{vtr > 200 kms~^) ~ 4%. This probability is consistent 
with both the currently observed MSPs and the local halo MSP number density nhainMRp ~ 
3 kpc ~^ relative to the disk MSP number density ndiskMSP ~ 40 kpc~^ (j Cordes fc Chernofi 
1993). 



In Figure [2], we consider the case that MSPs do not emit enough electrons and fit the 
7-ray haze data with the pulsed emission from MSPs and with DM annihilation products. 
In this case, we need a rather mild boost factor BF = 3. For the microwave haze data we 
need a significantly larger constant offset than in Figure [H and the shape of the microwave 
haze is not fit nicely due to the small production rate of e~^e~ by annihilating DM in this 
case. 

In both of the cases considered above the e'^e~ emission from MSPs and DM is too small 
to significantly contribute to local e~^e~ fiuxes. I n particular, we nee d to assume that the 



rising positron fraction in PAMELA experiment (lAdriani et al.l 120091) i s due to a di: 



astrophysical source, such as young radio pulsars (jHooper et al. 



2009 



Yiiksel et al. 



erent 



2009 



Profumo 



2008 



Malvshev et al.ll2009l ) 



In the case when e~^e emission from MSPs is small, one can also consider DM particles 



that predomi nantly annihilate into 



a DM model (lArkani-Hamed et al. 



epton s with large boost factors. As an exam ple, we study 



20091 1 used to explain the local e'^e fiuxes (jCholis et al. 



2009al ). The corresponding 7-rays signal is shown in Figure [31 The ICS photons have a fiat 
spectrum and cannot explain the bump around 4 GeV which can be produced by pulsed 
7-ray emission from about 20 000 MSPs in the Milky Way halo. Another possibility is to 
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Fermi Haze Template (|l| < 15, -30 < b < -10) 
Total 

XDM e*, ICS off ISRF 

XDMe*, FSRxlOO 

MSPs pulsed T, P,,=3.4x1 O^^erg s'^ 




Total 

XDM e*, Mjj=1.2 TeV, B.F.=100 
constant 
Haze Data : 



V = 23 GHz 



-t 



10 15 20 25 30 

Latitudinal Radial Distance from GC (degrees) 



35 



Fig. 3. — P ulsed 7-rav emission f rom MSPs and DM with e~^e annihilation mode. This DM model 
was used by 



Cholis et al 



(|2009al ) to fit the local e'^e" cosmic-ray experiments. The spectrum of 
ICS photons is too flat to reproduce the bump around 4 GeV. As a result, a contribution of 7-rays 
from a population of about 2 x 10^ MSPs in the Milky Way halo is necessary. 



consider a model with more than one species of DM particles (e.g., 



Cholis fc Weiner 



20091). 



In Figures H] and [5l we address the question whether the observationally unknown pa- 
rameters used in Figure [1] were fine tuned. In particular, we study the dependence on the 
index and cutoff of e^e~ injection spectrum from MSPs, on the scaling radius of the MSP 
distribution Rc, on the DM mass, and on the DM annihilation channel^. We show that 
rather large changes of these parameters do not change our conclusions. 

In this paper we have studied possible contributions to the 7-ray Fermi haze. We have 
shown that two major contributions may come from the millisecond pulsars and annihilating 
dark matter. The electron and positron emission from MSPs is uncertain and we considered 
two cases. In the first case, we assume that a large part of MSP spin-down luminosity is 
transformed in the e^e~ fiux. In this population of 3 x 10^ MSPs in the Milky Way 

halo can explain the 7-ray haze if about 50% of the MSPs spin-down is transformed into 
e~^e~ and about 10% is transformed in pulsed 7-rays. The synchrotron emission from MSP 



^ We show two cases of annihilation channels XX ~^ W^W~ and XX ~^ bb. The annihilation channel 
XX ^ T^T^ has a similar 7-ray spectrum to xX ~^ W^W^ but the number of high energy photons is lower 
by about a factor of 4. The channel xX ^ produces a softer 7-ray spectrum than xX ^ W^W~; as a 
result, a higher Mdm will be needed to account for the high energy part of the 7-ray haze. 
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electrons is consistent with the WMAP microwave haze data as welL Prompt 7-ray emission 
from a DM particle annihilating into W^W~ contributes significantly around 100 GeV. In 
the second case, we assumed that practically no MSP spin-down energy goes into e^e~ . The 
7-ray haze can be explained by the pulsed 7-rays from MSPs and a DM particle annihilating 
into W^W~ with a boost factor 3. However, the synchrotron emission (of e^e~ from our 
assumed DM model) is insufficient to explain the microwave haze data. Another alternative 
to explain the 7-haze in the case of small e^e~ emission from MSPs is DM particles with 
leptonic annihilation modes and boost factors of order 100 together with pulsed 7-rays from 
MSPs. The last scenario is also consistent with the WMAP microwave haze data. 

In general, the 7-ray haze, if confirmed, provides a strong evidence for new sources of 
7-rays additional to the standard ones, such as the 7r° production by cosmic rays and the ICS 
photons from the electrons accelerated by SNe in the Galactic plane. Possible new sources 
of 7-rays include DM annihilation and a large population of MSPs in the stellar halo of the 
Milky Way. The 7-ra y haze may als o suggest a recent AGN activity or an existence of a 



bipolar Galactic wind (jSu et al.ll2010l ). 
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Fermi Haze Template (|l| < 1 5, -30 < b < -1 0) 
XX-^W*W, Mj=500GeV, Pet=2.9x10%rg s 
XX^W*W", Mj=300GeV, Pet=2.7x10^''erg s 
XX^W+W", Mjj=200GeV, Pgt=2.5x10^''erg s 
xx^b6, M =300GeV. Pe==2.7x10'''erg s ' 




H Fermi Haze Tempiate (|i| < 1 5, -30 < b < -1 0) 

■ Ro=1 .0 l<pc, Pb==1 .8x1 O^^erg s"' , P.^5.8x1 O^^erg s"' 

■ Rc=1.5 kpc, PB==2.4x10^\rg s \ P.^8.2x1 O^^erg s"' 
Rc=2.0 kpc, PB==2.7x10^''erg s \ P.,=5.6x10^^erg s"' 
Rc=3.0 kpc, PB==3.5x10^''erg s"\ P.^1.2x10^Vg s"' 




E 
o 

> 
m 
O 

a- 



1 1 1 

: 1 — 1 — ^ Fermi HazeTempiate (|i| < 15, -30<b< -10) 


■ n=1.5, Pa±=2.6x10^''ergs"\P,^6.0x10^^ergs"' 

n=1.3, Pe±=2.7x10^Vg P^5.6x10^^erg s"' 

n=1.1, Pgi=3.3x10^^erg P,,=4.7x10^^erg s"' 














MSPs contNb.'^-N^ ^\ 




1 1 1 


1: 



S 1e-06 



1 1 1 

: 1 — 1 — ^ Fermi HazeTempiate (|i| < 15, -30<b< -10) 


■ E„„,_=100GeV, Pa±=1.8x10^''ergs"\P^8.0x10^^ergs"' ■ 

E„^,'=300GeV, Pa±=2.7x10^''ergs"\p,^5.6x10^^ergs"' " 


E„^,°=500GeV, Pgi=3.5x10^''ergs"\P,=3.6x10^^ergs"' ■ 










MSPs con1}i(b?'''~-~^^^^ 'Ov 




1 1 i\ 





Energy (GeV) Energy (GeV) 

Fig. 4. — Fits to 7-ray haze data from Milky Way halo MSPs and DM with varying parameters. In 
each plot we present both the combined contribution from MSPs and DM and the contribution of 
the component whose properties we vary in each case. The reference model is as shown in Figured) 
DM has a mass of 300 GeV and annihilates via XX. ^ W^W~ with (av) = 3 x 10^^^ cm^ s~^. The 
MSPs have distribution profile pmsp ~ + Rc)^, where Rc = 2kpc. The MSP e~^e~ reference 
injection parameters are rig = 1.3 and £^cutc = 300 GeV. The pulsed 7-ray emission parameters 



are n-, 



1.3 and E, 



cut^ 



XX 



W+W~ and 



-- 4 GeV. Upper left: varying DM mass Mdm = 200, 300, 500 GeV for 
300 GeV for xX ~^ Upper right: varying scaling radius Rc of MSPs 



distribution profile. Lower left: varying MSPs electron injection (and 7-ray) index (n = 
Lower right: varying MSPs electron injection cutoff E^q. On each plot, there are two lines of the 
same type: two solid lines, two dashed lines, two dotted lines, etc. The lower one has only MSPs 
or only DM contribution, while the upper one has MSPs plus DM contributions. 
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XX->W*W", Mj(=500GeV, Pg±=2.9x10^°erg s"' — 




-i 


XX->W+W", Mjj=300GeV, P3±=2.7x10^°erg — 




XX->W+W", Mjj=200GeV, P3±=2.5x10^°erg s"' ■■■ 






H^bb, Mjj=300GeV, P3±=2.7x10^°erg - - 






Haze Data ^ h 






V = 23 GHz 




± 
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Latitudinal Radiai Distance from GC (degrees) 



_ f 


n=1.5, Pe±=2.6x10^°ergs' 

n=1.3, P5±=2.7x10^°ergs' 

n=1.1,Pg±=3.3x10^*ergs"' - 

Haze Data r + i 


V = 23 GHz 

± 




\ ± 









10 15 20 25 30 

Latitudinal Radial Distance from GC (degrees) 



Rc=1 .0 kpc, P5I 
R^=1 .5 kpc, P5I 
R,=2.0 kpc, P,i 
R|.=3.0 kpc, P5I 



.2.4x1 O^'erg s ' 
.2.7x1 O^'ergs' 
.3.5x1 O^'ergs' 
Haze Data 




10 15 20 25 30 

Latitudinal Radiai Distance from GO (degrees) 



E,„,^=100GeV,P^i 
Ecui =300GeV, P^i 



=1.8x10^°ergs' 
=2.7x1 O^Vgs" 



Haze Data 



V = 23 GHz 



10 15 20 25 30 

Latitudinal Radial Distance from GO (degrees) 



Fig. 5. — Fits to TFMj4P microwave haze data from Milky Way halo MSPs and DM with the same 
parameters as in Figure [H 
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